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Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation is 
generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only 
for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields are 
found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian 
domains is made. 
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1 Introduction 

The large-scale magnetic fields in the Sun and other stars are 
thought to arise from a complex interplay between nonuni- 
form rotation, stratification, and convective turbulence (e.g. 
Ossendrijver 2003; Brandenburg & Subramanian 2005). 
According to mean-field theory of turbulent dynamos, one 
of the main mechanisms responsible for the generation of 
large-scale fields is the a effect, which is related to the ki- 
netic helicity of the flow in idealised cases (Moffatt 1978; 
Krause & Riidler 1980; Rudiger & Hollerbach 2004). In 
the simplest settings, kinetic helicity is produced by ro- 
tating inhomogeneous turbulence. The inhomogeneity can 
arise from density stratification, or, in the absence of strati- 
fication, from impenetrable boundaries. However, although 
many early numerical simulations of rotating convection did 
work as dynamos (e.g. Nordlund et al. 1992; Brandenburg et 
al. 1996), the magnetic field generated in these models was 
found to be mostly on small scales, with the typical spatial 
scale of the field of the same order as the typical size of 
convection cells. This posed a puzzle, given that the simu- 
lations included all the major ingredients (stratification and 
rotation) known to be necessary for the a effect to occur. 
Crucial insight into this problem was gained when the test- 
field method (Schrinner et al. 2005, 2007) was used to com- 
pute the relevant turbulent transport coefficients (Kapyla et 
al. 2009b). It was found that the a effect resulting from 
the relatively slow rotation used in the earlier studies was 
not enough to overcome turbulent diffusion. On the other 
hand, for more rapid rotation, turbulent diffusion was found 
to decrease which, according to mean-field models, would 
imply that large-scale dynamo action should become pos- 
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sible. This was recently confirmed by direct simulations in 
the same parameter regime (Kapyla et al. 2009a). Similar 
results had been obtained earlier in the geodynamo con- 
text (e.g. Jones & Roberts 2000; Rotvig & Jones 2002), 
but subsequent simulations by Cattaneo & Hughes (2006) 
suggested that large-scale field generation becomes impos- 
sible at larger magnetic Reynolds numbers. However, this is 
now believed to be a consequence of rotation being too slow 
(Kapyla et al. 2009a). 

Furthermore, during the last few years simulations have 
verified that shear can be the sole driver of large-scale 
dynamo action in non-helical turbulence (Yousef et al. 
2008a,b; Brandenburg et al. 2008). The shear-current effect 
(Rogachevskii & Kleeorin 2003, 2004) has been proposed 
as a mechanism responsible for the dynamo, but the details 
of the numerical results disagree in various respects with the 
theoretical expectations. This numerical finding provided a 
motivation to add large-scale shear into a convective system 
(Kapyla et al. 2008, 2009c; Hughes & Proctor 2009). As 
a result it was found that these models indeed yield strong 
large-scale dynamos, but their theoretical interpretation in 
terms of mean-field theory is not as clear as in the case 
where only rotation is present (Kapyla et al. 2009b). 

Local simulations can be used to demonstrate the ex- 
istence of a large-scale dynamo and the processes that are 
responsible for the generation of the magnetic fields can be 
more readily studied. However, they cannot be used to study 
the global behaviour of solar and stellar dynamos, which 
occur in spherical geometry and generate a much richer va- 
riety of physical effects. Furthermore, the large-scale flows 
that need to be imposed in the local models can arise self- 
consistently in the spherical models. 
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Simulations of convective dynamos in spherical geom- 
etry have a long history starting from the seminal papers of 
Oilman (1983) and Glatzmaier (1985), to the present (e.g. 
Brun et al. 2004; Browning et al. 2006; Brown et al. 2007, 
2009). These models are now quite advanced and, with suit- 
able boundary conditions, are capable of reproducing the 
solar internal rotation profile (e.g. Miesch et al. 2006). How- 
ever, no cyclic large-scale dynamos have so far been dis- 
covered in models of the solar convection zone (e.g. Brun 
et al. 2004; Browning et al. 2006). On the other hand, mod- 
els of rapidly rotating stars do exhibit significant large-scale 
magnetic fields (Brown et al. 2007, 2009), although none of 
these studies has yet been able to produce poleward prop- 
agating dynamo waves first observed by Oilman (1983), or 
anything akin to the solar cycle. 

Modelling a full convection zone in a spherical shell 
though clearly desirable, has an important drawback, with 
the spatial resolution achieved in such simulations being in- 
evitably lower than in models that encompass only a part 
of the shell. One can also argue that it may not be neces- 
sary to model the full latitudinal and longitudinal extents of 
the convection zone in order to capture the essential physics 
of the solar dynamo. Main arguments to justify this ap- 
proach are that the solar magnetic activity appears to be 
concentrated in a latitude range of ±40° around the equator 
and that the large-scale field is mostly axisymmetric. This 
has recently motivated simulations of forced turbulence in 
parts ('wedges') of full spherical shells (Brandenburg et 
al. 2007; Mitra et al. 2009a,b). These simulations aim to 
strike a reasonable compromise between the requirements 
of higher spatial resolution on the one hand and globality 
on the other. Convection simulations in wedge-shaped ge- 
ometry have already been performed by Robinson & Chan 
(2001) and DeRosa & Hurlburt (2003). Those simulations, 
however, did not include magnetic fields. In the present pa- 
per we report preliminary results of convection-driven dy- 
namos in spherical wedge-geometry. 

2 Model 

The simulations reported here are performed in spherical 
polar coordinates where (r, 9, cj)) denote radius, colatitude, 
and the azimuth, respectively (in some of our space-time di- 
agrams, e.g. Fig. |2l we use the latitude O = 90° — 9.) As 
our domain we consider, instead of a full spherical shell, 
a spherical wedge that spans 0.6R < r < R in radius 
(where R is the outer radius of the sphere which is also 
used as our unit length), — A8 < < AG in latitude, and 
< (p < A(j) in longitude. We then solve the following set 
of compressible magnetohydrodynamic equations for the 
magnetic vector potential A (which is related to the mag- 
netic field B by B = V X A), the logarithmic density In p, 
the velocity U, and the specific entropy s: 
DA 

— =UxB- fiovJ, (1) 
D\np 



— =g-2fnxU + -(Jx B-Vp+V ■ 2vpS) , (3) 
Dt p 



T— = -V- KVT + 2iyS' + ^ J^ 

Dt p p 



(4) 



where D/Dt = d / dt + U ■ 'V is the advective derivative, 
and J = Pq^'V X B is the current density. The vacuum per- 
meability is given by po, whereas magnetic diffusivity and 
kinematic viscosity are given by rj and i', respectively. The 
rotation vector is given by rj = r2o(cos 6', — sin6', 0), where 
rio is the constant rotation rate of the frame. For simplicity, 
we have ignored the centrifugal force, as it is a second or- 
der effect in fl; see Dobler et al. (2006) for a more detailed 
justification. The gravitational acceleration, g, is given by 

9 = — (5) 

where G is the universal constant of gravitation, M is the 
mass of the star, and r is a unit vector in the radial direction. 
The fluid obeys an ideal gas law p = pe{'y — 1), where p 
and e are the pressure and internal energy, respectively, and 
7 = Op I Cm = 5/3 is the ratio of specific heats at constant 
pressure and volume. The heat conductivity is given by K 
and the internal energy is related to the temperature via e = 
cyT. The specific entropy is given by 



Cy ln(p/po) - Cp ln(p/po) 



(6) 



where po and po are constants defined below. The tempera- 
ture can now be written as 



T = Tq exp 



7— + (7 
Cp 



1 



ln(^ 

Po 



where Tq = Po/[cvPo(7 - !)]■ 

The rate of strain tensor S is given by 



u 



(7) 



(8) 



where ';' denotes covariant differentiation (see Mitra et al. 
2009a, for more details). We consider a setup where the 
stratification in the region below r = 0.7i? is convec- 
tively stable, whereas the region above is unstable. The po- 
sitions of the bottom of the domain, bottom of the convec- 
tively unstable layer, and the top of the domain are given by 
(ri, r2, r^) = (0.6, 0.7, l)R. The setup yields a hydrostatic 
equilibrium temperature 



T(r) =T3 + 



cv(7 - l)im + 1) 



dr. 



(9) 



where T3 is the temperature at and where m — m{r) 
is chosen to have a suitable radial profile. This procedure 
gives a temperature stratification for which the logarithmic 
temperature gradient (not to be confused with the gradient 
operator V) is given by 



V = dlnT/dlnp = (to + 1) 



(10) 



Dt 



= -V-U, 



(2) 



The stratification is convectively unstable when V — Vad > 
0, where Vad = 1 — 1/7- Thus m can be considered as 
a polytropic index, with values to < 1.5 yielding an un- 
stable stratification. We choose to = 1 in the convectively 
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unstable layer (r > r2) and m = 3 in the stable layer be- 
low (ri < r < r2) with a smooth transition in between. The 
density stratification is obtained by requiring the hydrostatic 
equilibrium condition to be satisfied. 

The thermal conductivity is obtained by requiring a con- 
stant luminosity L throughout the domain via 



K 



iirr^dT/dr 



(11) 



In order to speed up thermal relaxation, we use m — 1 only 
for the thermal conductivity profile, whereas the actual ther- 
mal stratification of the thermodynamic variables followed a 
shallower profile with p ^ T^^. The initial non-convecting 
stratification is shown in Fig.[T] A fixed temperature gradi- 
ent is imposed at the base of the domain which leads to a 
constant flux of energy into the domain. At the outer radial 
boundary the temperature is kept constant and on the lati- 
tudinal boundaries we set dgT = so that there is no heat 
flux in or out of the domain. 

For the velocity, the upper radial and the latitudinal 
boundaries are taken to be impenetrable and stress free, 
whereas at the lower boundary rigid rotation is enforced by 
taking 



Ur 



C/0 = O, 
dUe 



0. 



dr 



dUg ^ Ue 
dr r 
_ Ug_ dU^ 
r ' dr 



^ (r = r,) 
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(12) 
(13) 
(14) 



Ur - 

dU 

where 9i = 90° - AG and 02 = 90° + AO. We use per- 
fect conductor boundary conditions for the magnetic field 
at the lower radial and latitudinal boundaries, and a normal- 
field condition on the radial top boundary. In terms of the 
magnetic vector potential these translate to 
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(r = ri), (15) 
(r = rg), (16) 
= 9i,92). (17) 



Finally, all quantities are taken to be periodic in the az- 
imuthal direction. 

The simulations were performed using the PENCIL 
COD^H which uses sixth-order explicit finite differences in 
space and a third-order accurate time stepping method (see 
Mitra et al. (2009a) for further information regarding the 
adaptation of the PENCIL Code to spherical coordinates). 
We use grid sizes of 64 x 192 x 128 (Runs A1-A4) and 
64 X 256^ (Run Bl) where the angular extent, 2 AO x Acf), 
is either 2 x 67.5° x 90° or 2 x 60° x 120°; see Table[T]be- 
low. Typical snapshots of the radial velocity Ur for the two 
domain sizes are shown in Fig. |2] 



T, p, p, 6 




Fig. 1 Initial stratification of temperature (solid line), 
density (dashed), pressure (dot-dashed), and the superadi- 
abatic temperature gradient (dash-triple-dotted). The sub- 
scripts b refer to the values at r = ri . 

2.1 Units, nondimensional quantities, and parameters 

Dimensionless quantities are obtained by setting 

R = GM^Po^Po = cp = fio = l, (18) 
where po is the density at ri . The units of length, time, ve- 
locity, density, entropy, and magnetic field are then given 
by 

[x]^R, [t] 



[U] = ^/GM/R 



^JWJgm , 

[p]=Po, W=cp, {B\ - ^poiioGM/R . (19) 

The simulations are now governed by the dimensionless 
Prandtl, Reynolds, Rayleigh and Coriolis numbers defined 
by 



Pr = 



Pm = 



7? 



Re 



Ra = 



Rm 



Co = 



Xo 

^rins 

GMjrs - r^Y 

vxoR^ 

2f7o 



PmRe, 



1 ds 
cp dr 



(20) 
(21) 

(22) 

(23) 



http : / /pencil -code . googlecode . com/ 



where xa = K/{pa^CY>) is the thermal diffusivity, fcf = 
27r/(r3 — is an estimate of the wavenumber of the 
energy-carrying eddies, p,n is the density in the middle of 
the unstable layer at rm — ^(^'s — ^2), and Wmis is the vol- 
ume averaged rms velocity. The entropy gradient, measured 
at rm in the non-convecting initial state, is given by 

L'}i) .Y^, ,24, 

cpdr/^^^ iip 

where Vm — (<91nT/91np)r,„, and Hp is the pressure 
scale height at r„i. 

The energy that is deposited into the domain from the 
base is controlled by the luminosity parameter 
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Fig. 2 (online colour at: www.an-journal.org) Radial velocity Ur near the upper boundary (r — 0.98i?) from 
Runs A4 (left) and Bl (right). The slices are projected onto a spherical surface which is tilted 15° towards the ob- 
server The dotted lines in the background show a grid in latitude and longitude in increments of 15°. See also 

http : / /www . helsinki . fi /'^kapyla /movies . html 



Table 1 Summary of the runs. Rm and Co are given in 
the saturated state of the dynamo. Pm = 1 in all runs. 



Run 


Rm 


Co 


Pr 


Ra 


£ 




Al 


36 





0.2 


1.2 ■ lO'' 


10"'' 


0.084 


A2 


71 


0.3 


0.2 


1.2 ■ 10^ 


10"^ 


0.084 


A3 


35 


1.3 


0.2 


1.2 ■ lO'^ 


10"^ 


0.084 


A4 


28 


2.5 


0.2 


1.2 ■ 10^ 


10"^ 


0.084 


Bl 


45 


5.5 


0.67 


2.6 ■ 10" 


1.2 ■ 10"" 


0.078 



Furthermore, the stratification is determined by the pressure 
scale height at the surface 

cv(7 - 1)73 



(26) 



GM/R ■ 

Similar parameter definitions were used earlier by Dobler et 
al. (2006). The parameter values used in our runs are listed 
in Table[T] 

3 Results 

3.1 Hydrodynamic solutions 

Before discussing our results with magnetic fields we briefly 
describe a limited range of hydrodynamic solutions. One of 
the most important reasons for using spherical geometry is 
that large-scale shear and meridional flows can be generated 
self-consistently when the overall rotation of the star is in- 
cluded. In such a setting the kinetic helicity is also likely to 
vary as a function of position. At large magnetic Reynolds 
numbers, this information can then be used as a proxy of 
the trace of the a tensor of mean-field dynamo theory (e.g. 
Radler 1980). In order to study the effects of rotation on the 
flow, we study here a range of models where the rotation 
rate is gradually increased. 



Figure [3] shows cross-sections of longitudinal and tem- 
poral averages of the meridional velocities, rotation profile, 
and kinetic helicity for the kinematic stages of the simula- 
tions, from three runs with varying degrees of rotational in- 
fluence. In the absence of rotation (Run Al), the meridional 
velocity components show a large-scale structure which is 
related to long-lived convective cells. However, as there is 
no mean rotation, no differential rotation is generated in this 
case, and, although the kinetic helicity oj ■ u shows fairly 
large peak values, there is no systematic component present. 
These results are in accordance with theoretical considera- 
tions for nonrotating convection. 

When rotation is relatively weak, with Co ~ 0.3 
(Run A2), very strong differential rotation and meridional 
flow is generated. The resulting differential rotation is anti- 
solar, i.e. the equator is rotating slower than the high- 
latitude regions. The difference between minimum and 
maximum rotation rates is remarkably large - about twice 
the rotation rate of the star. The meridional circulation is 
also very strong and oriented anti-clockwise. A similar re- 
sult has been reported by Robinson & Chan (2001) for their 
convection simulations in a wedge-shaped domain. How- 
ever, in their case the rotation profile changed to a more 
solar-like configuration when the longitudinal extent of the 
domain was doubled. We have not experimented with in- 
creasing the 0-extent, but we note that there are some other 
differences; for example, our domain is wider in longitude 
than that of Robinson & Chan (2001). The kinetic helicity in 
this case shows signs of large-scale structure: in the north- 
ern hemisphere there appears to be positive helicity near the 
equator and negative helicity at higher latitudes. A similar 
profile, but with opposite sign, is discernible in the south- 
em hemisphere. The strongest features appear near the lat- 
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Fig. 3 (online colour at: www.an-journal.org) From left to right: azimuthally averaged Ur, Ue,^ = / (r sin 9) 

and kinetic helicity lj ■ u from runs with Co = (top row, Run Al), Co « 0.3 (middle. Run A2) and Co « 2.5 (bottom. 

Run A4) from the kinematic stages of the simulations. Here Wrms is the rms-value of the vorticity, a; = V x u. 



itudinal boundaries and are possibly due to the converging 
large-scale flows that occur there. 

In the rapidly rotating case with Co « 2.5 (Run A4) 
the relative differential rotation is significantly reduced and 
equatorial acceleration is observed. However, the absolute 
value of dil/dr, relevant for the large-scale dynamo, is of 
the same order of magnitude, or somewhat larger near the 
equator, than in Run A2 with Co « 0.3. In agreement with 
earlier simulations of convection in rotating spherical shells 
(e.g., Miesch et al. 2000) the rapid rotation regime is dom- 
inated by large prograde 'banana cells' in the radial veloc- 
ity; see Fig. 121 In our most rapidly rotating models we also 
see nonaxisymmetric structures in the equatorial regions. 
Brown et al. (2008) described such structures as 'nests' of 



convection. The strength of convection in those structures 
varies periodically in time and the maxima coincide with en- 
hanced differential rotation similarly as reported by Busse 
(2002). The left panel of Fig. |2] shows indications of lo- 
calised convection whereas in the higher Rayleigh number 
model (Run Bl, right panel) these features have not yet de- 
veloped. We also find that these localised convection struc- 
tures are suppressed when magnetic fields of equipartition 
strength are present. The contours of constant il are aligned 
with the rotation vector, in accordance with the Taylor- 
Proudman balance 



dz 



0, 



(27) 
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where the z derivative is evaluated as 

d „ d sin 6' (9 

— =cose- — . (28) 

oz or r oO 

Only weak meridional circulation is observed and the flow 
consists of several small cells. This can also be attributed 
to the Taylor-Proudman balance and has also been found 
in mean-field calculations (Brandenburg et al. 1991). In the 
present models the Taylor-Proudman balance always dom- 
inates the rotation profile in the rapid rotation regime. Pre- 
liminary results from a run with roughly 20 times higher 
Rayleigh number (Run Bl) with Co « 5.5 gives qual- 
itatively very similar results. Reproducing the solar rota- 
tion profile self-consistently without resorting to an im- 
posed temperature difference between the equator and the 
pole (e.g. Miesch et al. 2006) or including unresolved turbu- 
lent convective heat fluxes (e.g. Durney & Roxburgh 1971; 
Kitchatinov et al. 1994; Rudiger et al. 2005) is likely to be 
difficult. 

We enforce rigid rotation at the lower boundary at r = 
0.6 R, although this is unphysical and not even enough to 
avoid spreading the differential rotation into the stably strat- 
ified layer below the convection zone between r = 0.6 R 
and r = 0.7 R (see also Miesch et al. 2000). This spreading 
is due to the viscosity which is much larger in our simula- 
tions than in the solar convection zone. Thus, no tachocline 
forms in our models. It is also likely that a self-consistent 
formation of a tachocline requires such high numerical res- 
olution that the problem cannot be tackled successfully in 
the foreseeable future. A more practical solution is proba- 
bly the damping of the velocities in the overshoot region, 
similarly as in Browning et al. (2006). 

We also find that in the most rapidly rotating case con- 
sidered here, the kinetic helicity has a well-defined pro- 
file with a mostly negative (positive) value in the northern 
(southern) hemisphere in the outer layers of the convec- 
tively unstable layer Furthermore, there is a sign change 
as a function of radius near the base of the convection zone. 
These results are in accordance with theoretical considera- 
tions and earlier results from local and global simulations 
(e.g. Miesch et al. 2000; Brun et al. 2004; Kapyla et al. 
2006). The simple first-order smoothing relation that is valid 
for isotropic turbulence states that a oc —a; • u which sug- 
gests that the a effect should in the present case have a pos- 
itive (negative) value in the northern (southern) hemisphere. 
This will be examined in more detail in the next section. 

3.2 Dynamo solutions 

We begin our study of dynamo excitation with a model with 
Rm « 36 that does not act as a dynamo in the absence of 
rotation (Co — 0). For Co w 0.2 the magnetic field shows 
exponential growth. The time averaged growth rate of the 
magnetic field, 

A= (dlnB,„,,/di)t, (29) 

increases with rotation for small values of Co, see Fig. 31 
although it shows a dip at Co « 2.5. 



0.020 




0.000 



Fig. 4 Growth rate of the magnetic field as a function of 
rotation. 
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Fig. 5 Root mean square value of the magnetic field as a 
function of time for runs with Co = (solid line. Run Al), 
Co w 0.3 (dashed. Run A2), and Co « 2.5 (dash-dotted. 
Run A4). 



In the slowly rotating cases (Co ^ 0.3) the saturation 
level of the magnetic field is rather low, Brms/Bcq « 0.1, 
where Bcq — \J Hopu'^^^^ is the equipartition field strength 
(see Fig.|5]l. Furthermore, there is a large-scale axisymmet- 
ric component to the field which is thus well described by a 
longitudinal average. In Run A2 with Co « 0.3 the large- 
scale field is approximately antisymmetric with respect to 
the equator and shows irregular sign changes on a long time 
scale (see the upper panel of Fig.|6]l. 

As the rotation rate is increased the qualitative charac- 
ter of the differential rotation changes from anti-solar to 
solar-like, as was seen in the previous section. A qualitative 
change occurs also in the behaviour of the dynamo so that 
for Co 1.3 (Run A3) magnetic fields of the order of the 
equipartition strength are obtained. Furthermore, the fields 
show a large-scale pattern that stays roughly independent 
of time, i.e. no sign changes are found. However, near the 
surface there are indications of a time-variable large-scale 
structure (lower panel of Fig. |6l). When the rotation rate is 
doubled (Co w 2.5, Run A4), the magnetic field indeed 
shows oscillatory behaviour near the equator, see Figs.Qand 
|9] with branches of magnetic fields appearing to move from 
the equator towards higher latitudes. We interpret the phe- 
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Fig. 8 (online colour at: www.an-journal.org) Azimuthally averaged magnetic field component i?^ from four times 
separated by At ~ corresponding to AtUims^f ~ 41, from Run A4 with Co « 2.5. 
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Fig. 6 (online colour at: www.an-journal.org) Az- 
imuthally averaged magnetic field component as a func- 
tion of time and latitude near the top (r = 0.98i?) of the 
convectively unstable layer for runs with Co w 0.3 (upper 
panel. Run A2) and Co w 1.3 (lower panel. Run A3). The 
horizontal white dotted line at = indicates the position 
of the equator. 



nomenon as a dynamo wave propagating towards the poles. 
More precisely, assuming that the sign of the a effect is op- 
posite to that of the kinetic helicity, we would thus expect a 
to have a positive (negative) value in the northern (southern) 
hemisphere. Furthermore, there is a positive dil/dr near 
the equator, which leads to a dft /dr > which, based on 
mean-field theory, should lead to poleward migration (e.g. 
Yoshimura 1975). We note that similar results were reported 
by Oilman (1983). However, a more detailed interpretation 
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Fig. 7 (online colour at: www.an-journal.org) Magnetic 
field components Br (top panel), Bg (middle), and (bot- 
tom) averaged over the azimuthal direction as a function of 
time and latitude near the top (r — 0.98i?) of the convec- 
tively unstable layer from Run A4. 

of the origin of the dynamo requires a better understanding 
of the turbulent transport coefficients. 

We note that in the kinematic stage of the simulation 
the period of the oscillations is well-defined and the sign 
changes occur nearly similarly on both hemispheres. In 
the saturated stage the oscillation period is expected to be- 
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Fig. 10 (online colour at: www.an-journal.org) Toroidal magnetic field i?^ near the upjjer boundary (r = 0.98i?) from 
the Runs A4 (left) and Bl (right). 
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Fig. 9 (online colour at: www.an-journal.org) Magnetic 
field component B^,, averaged over the azimuthal direction 
as a function of time and latitude from the bottom (top 
panel, r = 0.7R), middle (middle panel, r — Q.85R) and 
near the top (bottom panel, r = 0.98i?) of the convectively 
unstable layer from Run A4. 



come slightly longer when the implied turbulent diffusion is 
quenched by the magnetic field (e.g. Kapyla & Brandenburg 
2009). Such behaviour is indeed observed, but it appears 
that the two hemispheres begin to show a distinct asym- 
metry with different oscillation periods and magnitudes of 
the fields in the saturated state (see Figs.|7]|8] andfTOll. We 
note that hemispherical dynamos, where one of the hemi- 
spheres is almost devoid of magnetic fields, have been re- 



ported from rapidly rotating spherical shells (Grote & Busse 
2000; Busse 2002). These dynamos also exhibit poleward 
migration of activity belts similar to the results reported 
here. Furthermore, Fig. |9] shows that in the deeper layers 
the magnetic field has a very different structure with strong 
large-scale fields reaching much higher latitudes that can 
also change sign, but more irregularly and on a time scale 
that is much longer than that shown by the field near the sur- 
face. Since the differential rotation in these models is appre- 
ciable only near the equator, it is conceivable that there is a 
competition between different dynamo modes in the deeper 
layers and those at the higher latitudes and that the UTegular 
nature of the oscillations in the saturated state is caused by 
the interaction between these regions. 

4 Conclusions 

In this paper we have presented preliminary results of ro- 
tating turbulent convection in spherical wedge geometry for 
different rotation rates. We find that for small values of the 
Coriolis number the differential rotation is anti-solar, with 
the equator rotating slower in comparison to higher lati- 
tudes. When the Coriolis number is increased above unity, 
the rotation profile changes to a more solar-like configu- 
ration with faster rotation at the equator and slower rota- 
tion at higher latitudes. However, with more rapid rotation 
(Co 1) the Taylor-Proudman balance dictates the rotation 
profile with contours that tend to be constant along cylin- 
ders. 

The same setup is used to study the corresponding dy- 
namo solutions. The magnetic Reynolds number is chosen 
such that there is no dynamo in the absence of rotation. 
A dynamo is excited for sufficiently rapid rotation. Near 
the marginally excited state the magnetic energy is only 
roughly ten per cent of the kinetic energy, but a large-scale 
structure can already be observed. With more rapid rotation 
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(Co ^ 1), strong large-scale magnetic fields are obtained. 
In this case the field is non-oscillatory. Increasing the rota- 
tion rate further changes the solution such that an oscilla- 
tory magnetic field is observed near the equator in the upper 
parts of the convectively unstable layer The magnetic fea- 
tures appear at low latitudes and propagate poleward, i.e. 
opposite to the behaviour seen in the Sun. However, the ki- 
netic helicity in this run suggests that the q effect is positive 
(negative) in the northern (southern) hemisphere. This to- 
gether with the positive radial gradient of near the equator 
is consistent with a dynamo wave propagating poleward. 

Our model is in some ways similar to that of Gilman 
(1983) which also omitted the regions near the pole and 
used the same boundary conditions for the magnetic field. 
The main difference between our models is that we include 
density stratification and add a lower overshoot layer Many 
of the results, e.g. a poleward propagating dynamo wave in 
our most rapidly rotating runs, are qualitatively similar to 
those found by Gilman (1983). However, we obtain a deep- 
seated non-oscillatory field and see large differences in the 
activity levels on different hemispheres which were not re- 
ported by Gilman. The differences between the setups could 
explain some of the discrepancies of the results. 

Our ultimate aim with this model is to study dynamo 
processes that operate in solar and stellar interiors. In or- 
der to achieve this goal we need to find a way to break 
the Taylor-Proudman balance. A possible solution would 
involve the inclusion of subgrid scale anisotropic turbulent 
heat transport in the model. It is also unclear how the in- 
troduction of latitudinal boundaries affects the large-scale 
behaviour of the solutions. Issues related to these problems 
will be dealt with in future pubUcations. 
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